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ABSTRACT 

The binary confusion noise spectrum at LISA band depends strongly on observational 
period and abundance of Galactic close white dwarf binaries (CWDBs). We have 
investigated how the number of the resolved Galactic CWDBs varies with operation 
period of LISA, and found that the resolved number would typically grows by a factor 
of 5 when the operation period increases from lyr to lOyr. We have also made a similar 
estimation for number of CWDBs whose chirp signal can be measured in matched 
filtering analysis. 

■ Key words: gravitational waves-binaries: close-white dwarfs. 

m 
<N 
O . 

. l INTRODUCTION 



o 



o 



If things go well, the Laser Interferometer Space Antenna (LISAQ) will be launched within ten years. LISA and ground- 
based detectors (such as, TAMA300, GEO600, LIGO and VIRGO) will bring us fruitful information of our universe and era of 
gravitational-wave astronomy will really start. For example, using LISA, we might detect merging massive black holes (MBHs) 
with significant signal to noise ratio (SNR), or in-spiraling compact stars (such as, white dwarfs, neutron stars, stellar mass 
black holes) around MBHs. These are very exciting phenomena and we might confirm existence of black holes, make stringent 
Q$ . tests of general relativity and measure various interesting parameters of MBHs (Bender et al. 1998). However, event rate of 
^ ' MBH-MBH binaries is highly unknown (e.g. Haehnelt 1994) and gravitational waves from in-spiraling compact stars around 
• i-H , MBHs might be too much complicated to be detected by usual matched filtering technique (Bender et al. 1998 and references 
' therein). 

Galactic binaries, such as, neutron star binaries, cataclysmic binaries, close white dwarf binaries (CWDBs), are guaranteed 
sources for LISA (Mironowskii 1965, Evans, Iben & Smarr 1987, Hils, Bender & Webbink 1990). The Galactic CWDBs are 
expected to be the dominant one in the frequency region from ~ 10 _3 Hz up to several 10 _2 Hz (Bender et al. 1998). At present 
only less than ten CWDBs have been optically detected, but LISA would find thousands of CWDBs with one year integration 
(Hils & Bender 1997, Bender et al. 1998). Observational analysis for abundance and spatial distribution of Galactic CWDBs 
would bring us important clues to understand formation of binary stars and structure of our galaxy (e.g. Ioka, Tanaka & 
Nakamura 2000, Hiscock et al. 2000, see also Oppenheimer et al. 2001). It should be also noted that CWDBs are regarded as 
one of the likely progenitors of type I supernova (Iben & Tutukov 1984, Branch et al. 1995). Thus gravitational waves from 
Galactic CWDBs, one of the guaranteed sources of LISA, have rich scientific contents. 

An interesting feature of Galactic CWDBs is that they would also become an serious background noise (confusion noise) 
for measurement of gravitational wave (Evans et al. 1987). Properties of its spectrum as well as the estimation errors for 
parameters of resolved binaries depend strongly on observational period. But operation period of LISA is not determined 
definitely at present. The Pre-Phase A Report (Bender et al. 1998) of LISA states "This (a 1-year-long observing period) is 
reasonable length of time, but not the maximum: the nominal mission lifetime is 2yr, but in principle it might last as long 
as a decade.". Thus it seems interesting to investigate ahead of time how scientific impacts on CWDBs would change with 
operation period of LISA. We mainly investigate the number of the resolved CWDBs for various observation periods. This 
quantity would be most fundamental for observational analyses described in the last paragraph. 



* http://lisa.jpl.nasa.gov 
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This article is organized as follows. In §2 we mention gravitational wave from CWDBs and a model of their spatial 
distribution in our galaxy. In §3 we determine the confusion noise curve as a function of the observational period and the 
abundance of Galactic CWDBs. We also discuss signal to noise ratio and the parameter estimation errors for binaries in 
matched filtering analysis. Then in §4 we evaluate number of resolved CWDBs for various set of parameters. §5 is devoted to 
a brief summary. 



2 GRAVITATIONAL WAVE FROM GALACTIC CWDBS 
2.1 Chirp Gravitational Wave 

Let us consider a binary (masses Mi and M2) with a circular orbit. The characteristic amplitude Ka of gravitational wave at 
a frequency / is evaluated by Newtonian quadruple formula as (Thorne 1987) 

/2 £5/3 



flA 



= 8 (S) 12 ^^ 2/3 w 

Mc \ 5/3 /ioo P cw / y/3 

O.3M J \ r ) VlO- 3 lW ' 1 ; 



= 1.2 x 10" 21 



where M c = (MiM 2 ) 3/5 (Mi + M 2 )" 1/5 is the chirp mass of the system and r is the distance to it. In the above equation 
we have taken angular average of wave amplitude with respect to orientation of sources (Thorne 1987). The time profile of a 
nearly monochromatic wave is expressed as follows 



h(t) = Ha cos 



2- ( ft+ f -t 2 



(3) 



Here <f>o is an integral constant. For a purely monochromatic wave we have / = 0, but energy loss due to gravitational radiation 
or other effects, such as, accelerating motion of the binary or mass exchange between binary stars, make modulation of the 
wave frequency / ^ 0. We denote / 7^ as follows 

/=(/W + (/Lw> w 

where the first term of the right-hand-side represents effect due to gravitational radiation. When the wave frequency / changes 
significantly within an observational period T, namely 

fT/f > 0(1), (5) 

we might separate these effects in matched filtering analysis using their frequency dependence. But this is not the case for 
our analysis for low-frequency waves. Thus we try to measure / as the sole parameter for the frequency modulation. Effects 
due to motion of LISA would be commented in §3.2. 

Our main target is Galactic CWDBs. A CWDB is expected to have an almost circular orbit as a result of spiral-in process 
during its common envelope phase, in contrast to a double neutron star binary (Ignatiev et al. 2001). Before the less massive 
star fills its Roche lobe, gravitational radiation reaction is the dominant process for evolution of orbital parameters of a binary 
(see e.g. Webbink 1984). The energy loss due to gravitational radiation is given by Newtonian quadruple formula and changes 
the wave frequency as follows (Thorne 1987) 

« OT - ^^/""M!" (6) 



= ^ '"(^)"'mr-* 



CWDBs are mainly divided into three segments (Webbink 1984, Iben & Tutukov 1984, Tutukov & Yungelson 1986). They 
are He-He systems with total mass 0.5 — 0.75M Q , He-CO systems with total mass 0.75 — 1.45M , and CO-CO systems with 
total mass 1.45 — 2.4M . Webbink (1984) discussed that these three systems have similar abundance (see also Branch et al. 
1995). In this article we do not distinguish them but study a single component with chirp mass M c = 0.3Mq. More detailed 
treatment would not significantly change our main conclusions. 



2.2 Number and Spatial Distributions of Galactic CWDBs 

The abundance of Galactic CWDBs has not been observationally clarified (Hils et al. 1990, Marsh 1995, Marsh, Dhillon, & 
Duck 1995, Bender & Hils 1997, Knox, Hawkins, & Hambly 1999). In the frequency region where they evolve by gravitational 
radiation reaction, the number of CWDBs per unit frequency dN/df becomes 
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dN fdfy 1 11/3 

from equation (Q). The above distribution can be also expressed in an integral form as N(> /) oc f~ s ^ 3 . 

We use estimation for the abundance of Galactic CWDBs given in Hils & Bender (2000) as a reference. This estimation 
is 10% of the theoretical estimation by Webbink (1984). Hils & Bender (2000) commented that the error in this estimation 
is believe to be no more than a factor of ten in either direction. The number distribution dN/df corresponding to this 10% 
estimation is a factor 10 lower than that given in their figure 4 which represented the full 100% estimate based on Webbink 
(1984). Thus we have an approximation for this 10% estimation as 

at the frequency region / <; 10~ 3 Hz in interest. In this article we use the following non-dimensional quantity A r 282 to 
characterize the abundance of Galactic CWDBs 

1 dN 



N 



(10) 



/ = 10 -2.82 Hz 



" (3.2 x Wmz- 1 ) df 
or we have 

dN / f \ ~ I 1 / 3 

^-3.2xl0 7 iV 282 ( IF 4 1 -) Hz-. (11) 

For example we have N2S2 = 1 for distribution (^). Webbink & Han (1998) discussed CWDBs with various models of binary 
evolution. Their "standard model" corresponds to N2S2 ~ 2.0. Nelemans, Yungelson, & Portegies Zwart (2001) obtained a 
similar results. Note that (3.2 x 10 7 ) _1 Hz is the frequency bin for lyr integration. Thus A^82 represents the mean number of 
CWDBs within lyr" 1 bin at / = 10 _2 - 82 Hz. 

Abundance of extra- Galactic CWDB also becomes important to discuss the binary confusion noise in the next subsection. 
We scale their abundance by the same parameter JV282 defined for the Galactic ones (see e.g. Bender & Hils 1997 for discussion). 

Next we briefly discuss the spatial distribution of Galactic CWDBs. We use the standard exponential disk model 

p{R,z) = poexp (~^r) ex P ' ( 12 ) 

where (R, z) is the Galactic cylindrical coordinate. We fix the radial scale length Ro — 3.5kpc and the disk scale height 
Zo = 90pc (Hils, Bender & Webbink 1990), and assume that the solar system exists at the position R = 8.5kpc and \z\ = 30pc. 
As the amplitude of the gravitational wave is inversely proportional to the distance r, distribution of the CWDBs' distances r 
from the solar system becomes important. With our model parameters, ~ 90% of the Galactic CWDBs are within r < 18kpc, 
and ~ 10% of them are within r < 5.0kpc. The number of binaries within distance r changes from r3 ( r: smaller than 
disk thickness ~ lOOpc) to r 2 (r: smaller than size of Galaxy ~ lOkpc) and r ° ( r: larger than ~ lOkpc). 



3 MATCHED FILTERING ANALYSIS 
3.1 Confusion Noise 

The signal of a detector s(t) = h(t) +n(t) contains both the true gravitational wave signal h(t) and the noise n(t). We assume 
the stationary noise and define its power spectrum S„(f) by 

(n(f)n(f')) = \s(f-f')S n (f), (13) 

where n(f) = f^ a e 2n '-' t n(t)dt is the Fourier transformation of the noise n(t). This definition of S n (f) corresponds to the 
one-sided spectral density (e.g. appendix A of Cutler and Flanagan 1994). For notational simplicities we have represented the 
frequency space in the continuum limit (see e.g. Schutz 1997). 

At the LISA band the noise spectrum S„(f) is constituted by two terms; (i) the instrumental (detector's) noise Sms(f) 
and (ii) the confusion noise S CO n(f) that is effectively caused by unresolved sources 

S n (f) = Sins{f) + Scon(f)- (14) 

The frequency bin 8f for an observational period T is simply given as 
T x 



8f = T = 3.1 x 10 ( y— J Hz. (15) 
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The number dN/df of Galactic CWDBs per unit frequency increases at lower frequency as shown in equation (g). Roughly 
speaking, when the frequency bin (|l^) is occupies by more than one Galactic CWDBs, the Galactic confusion noise S CO n,G(f) 
becomes important (e.g. Bender & Hils 1997). At higher frequencies there is no more than one Galactic binary and the 
confusion noise S CO n,eG(f) is determined by the extra-Galactic binaries. The shape of the confusion noises (including the 
position of the transition frequency ft where the Galactic noise becomes important) is determined by the observational period 
T and the number of binaries. We discuss them with a method similar to Hiscock et al. (2000) who studied dependence of 
the noise spectrum on the abundance and spatial distribution of CWDBs (see also Phinney 2001). 

Bender & Hils (2000 and references therein) discussed that the confusion noise in the region 10~ 4 < / < 10~ 15 Hz is 
dominated by CWDBs. They found that in some scenarios (Tutukov & Yungelson 1996, see also Iben & Tutukov 1991) the 
abundance dN/df of the Helium Cataclysmics (HeCVs) exceeds that of CWDBs, but amplitude of confusion noise is expected 
to be dominated by CWDBs because of their larger chirp masses than that of HeCVs. Thus in our present analysis we 
basically do not consider the effects of HeCVs and study the spectrum S CO n(f) as a function of T and abundance of CWDBs 
(characterized by N2S2). 

The extra-Galactic noise S CO n,eG(f)df has the following relation 

S CO n,eG(f)df oc (number of binaries in / ~ / + df) x (amplitude of individual source) 2 . (16) 
From equations (Q) and (^) we have 

SconMf) OC iV 282 /- 7/3 . (17) 

The Galactic confusion noise S CO n,G(f) at frequency smaller than the transition region around the frequency ft has the same 
functional shape S CO n,G(f) oc A r 282/ _7//3 . The ratio of two amplitudes yj Scon,c(f)/ S CO n,eG(f) is estimated to 3 ~ 10 and the 
uncertainty is mainly related to that of star formation rate at high redshift (Kosenko & Postnov 1998, Schneider et al. 2001). 
Next we evaluate the transition frequency ft as a function of the parameters -/V282 and T 



dN 
df 



Sf = Y, (18) 

ft 

where Y is a constant of order unity. From equation ( pd| ) we obtain 

ft oc (7V282T- 1 ) 3 / 11 . (19) 

Now we can make the confusion noise spectra S CO n(f, N282, T) for various parameters (N2S2,T) using a functional form 
S con (f, N282,o, To) given for a specific choice of parameters iV282,o and To. We connect two (Galactic and extra-Galactic) 
power-law functions oc A r 282/ _7/ ' 3 around the characteristic frequency ft and straightforwardly obtain the following relation 

S con (f,N2S2,T) = ( ( l^) 7/S S con (f ( f jA ,N ,T ) , (20) 



where we have 

ft,0 _ I N282,oT ~ 



-^282,0 / V ft J V V ft 



3/11 



(21) 



ft V ^282 T- 1 

from equation dlj| ) 

For the "base" function of the confusion noise S con (f,N2S2,o,To) we use the result of Hils & Bender (2000) given for 
parameters N2S2 = 1 and To = lyr (see their figure 5). The instrumental noise S n (f) is obtained from their same figure. In 
figure 1 we show the noise spectrum for various choice of parameters (N2S2,T). 

Bender & Hils (1997) gave confusion noise curves S CO n(f) due to CWDBs for two different A^282- We find that those 
two curves are reasonably reproduced by each other with the above procedure. Our method would be quantitatively valid at 
frequency 

/ rp \ -3/11 

/>3.0x 10- 3 (7V 2 82) 3/1 M — J =ft, (22) 

where the confusion noise is determined by extra-Galactic binaries. We have put Y = 0.1 in equation ([l8]) for definition of 
ft- The Galactic binaries with / > ft are not overlapped in the frequency bin T" 1 and, in principle, simple to analyzed. In 
contrast these with f < ft can be detected only when they are close to us and have significant wave amplitude above confusion 
noise. Their analysis would not be straightforward. Here we effectively discuss them using the confusion noise spectrum. 

In figure 1 the frequency ft is shown with filled circles. Note that the Galactic confusion noise at the transition region for 
T = lyr vanishes in the case of T = lOyr. Comparing T = lyr and T = lOyr we find that y S CO n(f) of the latter is smaller 
about a factor of ten than the former at / ~ 0. 0018Hz in the case of -/V282 = 1. Note that this factor is determined by the 
quantity S co „,g(/) / S CO n,eG(f) as discussed before. 
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Figure 1. The noise spectrum is presented in the form of */ S(f). The instrumental noise of LISA is presented with bold thick line 
(Hils & Bender 2000). Other lines are confusion noise with various set of parameters N2S2 and T (in units of year). The thin solid line 
corresponds to (N2$2,T) = (1, 1), thick sold line to (1,10), long-dashed line to (0.1,1) and short-dashed line to (10,1). The filled circles 
represent the frequency ft defined in eq. 



3.2 SNR and Estimation Errors for a Chirp Signal 

In this subsection we briefly discuss signal-to-noise ratio (SNR) and the parameter estimation errors for gravitational wave 
signal in matched filtering analysis (Cutler & Flanagan 1994). First, we define an inner product of two quasi-periodic waves 
g(t) and k(t) around a frequency fo as follows 

(9\k) = ^|tt / g(t)k(t)dt. (23) 
The signal to noise ratio of a gravitational wave h(t) is given as 

SNR = (h(t)\h(t)) 1/2 . (24) 

When the wave form h(t) is characterized by some parameters Ai, magnitude of the estimation errors AAi is given by the so 
called Fisher information matrix Vij as 



<AA,AA,)= -I- =r- (25) 



/ dh 


dh \ 


\d\i 


dxj 



Position and orientation of LISA change in time. This causes (i) frequency modulation due to Doppler effect and (ii) variation 
of detector's sensitivity due to its rotation (Bender et al. 1998). In reality we need to fit the direction of a source in the 
matched filtering analysis (Peterseim et al. 1997, Cutler 1998), but this is a very troublesome task. In the present analysis we 
use an angular averaged sensitivity (effectively a factor of %/B" degradation, Thorne 1987) and do not try to fit the direction. 
We evaluate the estimation errors for the three parameters of wave form Xi = {/, / , 0o} in equation (^). The chirp signal / is 
very important from astronomical point of views (see e.g. Schutz 1986, 1989). This parameter is related to secular effects of 
gravitational wave during the whole observational period T and distinct from annual effects related to the direction of a source 
described above. However, we should notice that due to correlation in the Fisher information matrix the actual parameter 
estimation errors AAi would be worse than our results. Thus we must be cautious to discuss the measurement error of the 
chirp signal /. 

Now we can calculate SNR and the estimation errors AAi for Ai = {/, /, 0o}- Using the time function (^|) and equations 
( poj ) and (Q) we obtain the following results for a binary 

SNR= —^=T 1/2 , (26) 
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and the information matrix Yij gives the root-mean-square values of the errors as 

A *» - ws- < 29 > 

For comparison with equation ([?]) we can write equation ( |27| ) in the following form 



A/ = 4.3xl0--f^" 1 fJL 
J V 100 / WOyr 



-2 



(30) 



At a frequency / where the noise S n (f) does not depend on the observation period T (as in the case of / S> ft- see figure 
1), the above quantities have the following relations 

(31) 
(32) 
(33) 
(34) 

Relation (J3TJ) is a well known result. Note that the error ( |32| ) for the chirp signal becomes significantly smaller for longer 
integration time T. Considering the correlation of the Fisher matrix discussed before, the amplitude of error Af might be 
larger than our estimation (p7j). But the asymptotic time dependence A/ oc J 1-5 / 2 would be same for integration time T 
much larger than time scale (lyr) of annual modulation due to rotation of the detector. 



SNR 


■oc 


T 1/2 , 


A/ 


oc 


ji-5/2 


A/ 


oc 


ji-3/2 


A</> 


oc 


ji-l/2 



4 RESOLVED CWDBS 

In figure 2 we show the number of the Galactic CWDBs resolved with SNR > 10. We take our model parameters for abundance 
of Galactic CWDBs at -/V282 = 0.1, 1.0 and 10 and for observational period at T = lyr and lOyr. We present the number of 
the resolved CWDBs within frequency bin fa < f < 1.26/o(= 10 1,/10 /o) and explain our results mainly using the figure (the 
upper right panel of Fig. 2) given for N2&2 = 1.0. 

Firstly, we should notice that all the Galactic CWDBs with / ^ 3 x 10~ 3 Hz are resolved with SNR > 10. Thus we 
have N(fo < f < 1.26/o) oc fodN/df\f oc fg 8 for CWDBs at higher frequency. Secondly, according to our analysis based 
on the confusion noise spectrum, the number of resolved CWDBs increases significantly at frequency region dominated by 
Galactic confusion noise (/ < 10 _3 Hz) for T = lOyr comparing with T = lyr. In figure 3 we plot the distances r of CWDBs 
measured with SNR = 10 in the case of N2S2 = 1. At / < ft the distance r suddenly becomes small due to the Galactic 
confusion noise. As the observational period increases, the frequency ft decreases and we can observe more CWDBs. In table 
1 we show the total number of resolved CWDBs for various parameters. We also count CWDBs only with / > ft. We can 
understand that resolved CWDBs mainly belong to this group that do not overlap with other Galactic CWDBs. The total 
number for T = lOyr is 4 ~ 6 times larger than for T = lyr. As expected, this improvement is more significant for larger iV282- 
Note that the peak of the distribution N(fo < f < l-2/o) corresponds to the characteristic frequency ft for N282 = 1 and 10. 
This means that we have N(> ft) oc / t ~ 8//3 ^ J ,8 / 11 ; and the factor ~ 5 mentioned above is roughly given by 10 8//n ~ 5. 

A He+He white dwarf binaries coalesce at / ~ 0.015Hz and the He+CO at / ~ 0.03Hz. The coalescence frequency of 
CO+CO white dwarf is close to 0.1Hz (Bender & Hils 1997). For simplicities we do not take into account of these cut-off 
frequencies. 

Cutler (1998) studied LISA's angular resolution AQ for binaries with monochromatic gravitational waves (see his table 
1). The resolution depends on the orientations of binaries as well as the directions to them. Typical estimation errors with 
lyr observation become Ail ~3x 10 _2 [sr] at / = 10" 3 Hz and Afi ~lx 10 _3 [sr] at / = 10~ 2 Hz (both for binaries with 
SNR =10). 

Next we investigate the number of CWDBs whose chirp signal due to gravitational radiation reaction (eq.Q) can be 
measured. We define a parameter C by the following equation 

(f)aw- C ( A -f)- (35) 
For reference we take C > 10 as a detection criteria for the chirp signal by the gravitational radiation reaction. This threshold 
means that the signal of magnitude (f)aw (eq. M) is measured within 10% accuracy in matched filtering analysis. The number 
distribution JV(/o < / < 1.26/ ) of CWDBs with C> 10 is shown in figure 2. Note that CWDBs with C > 10 has SNR > 10 
for reasonable choice of T. We can directly confirm this using equations (Q) and (^). In table 1 we show the total number 
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Table 1. number of the resolved Galactic CWDBs 



(iV282,T[yr]) 


(0.1,1) 


(0.1,10) 


(1.1) 


(1,10) 


(10,1) 


(10,10) 


SNR > 10 


350 


1460 


1950 


10700 


4810 


28500 


SNR > 10, / > ft 


325 


1420 


1610 


9340 


3810 


22500 


C> 10 


0.34 


13 


3.4 


130 


34 


1200 



-^282 represents abundance of the Galactic CWDBs as defined in equation (ftd) and T is the operation period of LISA in units of year. 
C is the resolution threshold for chirp signal (eq. |35|). 



N c hirp with C > 10. The number of CWDBs with C > 10 increase significantly for long term integration. The total number 
of CWDBs for T = lOyr is ~ 40 times larger than T = lyr. 

As explained earlier, the actual estimation error A/ might be larger than our evaluation ( p7| ) due to correlation with other 
errors, such as, angular position of the source. Here we discuss how the total number N c hi rp depends on the observational 
period T assuming the following relation from equation (Ji?)) 

A/ = Ax (SNR^T- 2 , (36) 

with a normalization constant A that would be larger than 3\/5 /n in equation (E^). This parameter A is effectively absorbed 
to the threshold value C. We can expect following two points from figure 2 (i) the total number N c hir P would be mainly 
determined by f max where the distribution fdN/df takes maximum value, and (ii) it also departs from the asymptotic 
behavior oc N282.f~ 8 ^ 3 at high frequency region. The former point (i) represents that number N c hi rp of resolved binary 
becomes 

N chlrp oc AT 282 /-^ 3 . (37) 

As the effects of the confusion noise is very small around / ~ fmax, the total number N cn i rp simply scales as N c hi rp ^ N2&2 
in contrast to the previous case of the SNR threshold. The latter point (ii) means that the distance r of CWDBs for a given 
threshold C corresponds to r ~(size of Galaxy) at / = fmax- Then we can relate T and f ma x using equations (^) and ( |35| ) 
as follows 

J.-5/2 

J-2/3 
J max 



\J Sn(fmax) OC fmax 1 (38) 



and obtain 

fmax — fmaxSni^fmax^) OC T , (39) 

where we have used S n (f) — Sins(f) 25 / m l ne frequency region relevant for the present analysis (see figure 1). Using the 
relation (^) we finally obtain 

Nchirp OC A^282T 20/13 . (40) 

Thus a long observational period T is crucial to increase the total number N cn i rp - 



5 SUMMARY 

The Galactic close white dwarf binaries (CWDBs) are one of the guaranteed sources of LISA. Gravitational wave astronomy 
for the Galactic CWDBs would bring us important observational facts to understand binary formation, galactic structure, 
progenitor of type I supernova and so on. But gravitational waves from the Galactic CWDBs would also become a serious 
noise (called confusion noise) below the typical frequency /t(oc T~ 3//11 ) (eq.|E3]) where the effective frequency bin T -1 is 
occupied by more than one Galactic binaries (T: observational period). Apparently the confusion noise spectrum depends 
strongly on the observational period T. 

We have investigated number of the resolved Galactic CWDBs as a function of T using the exponential disk model. We 
found that the number of the resolved CWDBs with SNR > 10 increases y 8 / 11 as a function of observational period 
and it becomes about a factor of 5 larger for T = lOyr comparing with T = lyr. We also studied the number of CWDBs 
whose chirp signal can be measured with matched filtering method. The chirp signal is one of the most fundamental signal for 
gravitational wave astronomy. Using a rough estimation we have shown that the number would grow strongly ( oc 7-20/13) 
with period T. From lyr to lOyr operation the number can increase a factor of 40 for a typical model of Galactic CWDBs. 
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Figure 2. Number distribution of resolved Galactic CWDBs within frequency bin /o < / < l-26/o for parameters -/V282 = 0.1, 1 and 10. 
Solid lines represent results for T = lyr and dotted lines for T = lOyr. Thick lines are number of resolved CWDBs with SNR > 10, and 
thin lines with R > 10 (chirp signal). 




0.0001 0.001 0.01 0.1 



Figure 3. Distance r of the Galactic CWDBs resolved with SNR > 10. We fix the abundance by iV282 = 1- The sold line represents 
results for T = lyr and the dotted line for T = lOyr. 
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